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We present the first search for supersymmetry (SUSY) in Zy final states with large missing 
transverse energy using data corresponding to an integrated luminosity of 6.2 fb _1 collected with 
the DO experiment in pp collisions at yfs — 1.96 TeV. This signature is predicted in gauge-mediated 
SUSY-breaking models, where the lightest neutralino x? is the next-to-lightest supersymmetric 
particle and is produced in pairs, possibly through decay from heavier supersymmetric particles. 
The x? can decay either to a Z boson or a photon and an associated gravitino that escapes detection. 
We exclude this model at the 95% C.L. for SUSY breaking scales of A < 87 TeV, corresponding to 
neutralino masses of M(x?) < 151 GeV. 
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Gauge-mediated SUSY breaking (GMSB) 1] is a well- 
motivated model for physics beyond the standard model 
(SM). In GMSB models, SM gauge interactions serve as 
the messengers of SUSY breaking and thereby the masses 
of the SUSY partners of SM particles are connected to the 
strength of their gauge interactions. Assuming R-parity 
conservation, SUSY particles are produced in pairs, each 
decaying to lighter states which always include the next- 
to-lightest supersymmetric particle (NLSP). The final su- 
persymmetric decay of the NLSP to SM particles and the 
nearly massless gravitino G provide the typical signature 
used in GMSB searches. A recently formulated model- 
independent framework for gauge mediation is discussed 
inRef. @. 

The CDF, DO, ATLAS, CMS, and HI Collaborations 
have all searched for GMSB neutralinos \i m the 7G+7G 
(and single 7G) final state, assuming that the y? is the 
NLSP and decays promptly producing a photon [3j-|5j . In 
this Letter, we present the first search for the ZG + 7G 
final state. The minimal GMSB model we consider is 



"Model Line E" of Ref. 16] which is characterized by six 
parameters: the effective SUSY-breaking scale A which 
is varied in the following, the number of sets of messen- 
ger particles which is set to 72,5 = 2, the ratio of the 
Higgs vacuum expectation value which is chosen to be 
tan/3 = 3, the mass of the messenger particles which is 
selected to be M/A = 3, the Higgs sector mixing pa- 
rameter [i which is taken as // = (3/4)Afi where M\ is 
the hypercharge gaugino mass, and the parameter C grav 
which is linearly related to the gravitino mass and is set 
to Gg ra v = 1 [7J. I n this model Xi decays with substan- 
tial branching fraction to ZG, as well as to 7G, thereby 
providing a promising experimental signature for the dis- 
covery of the X? NLSP in the ZG + 7G final state. The 
gravitinos escape detection, leading to a Z7 final state 
with large missing transverse energy, f, T . We report a 
search for these events in pp collisions recorded with the 
DO detector [8(] at the Fermilab Tevatron Collider. 

The final state for this analysis contains a Z boson 
decaying to e + e~ or /i + jjT , a photon of large transverse 



energy, and large f T . The data have been collected using 
a set of inclusive electron or muon triggers, corresponding 
to an integrated luminosity of 6.2 ± 0.4 fb _1 [9]. The 
triggers have about 100% (78%) efficiency for signal in 
the ee7 (/J./J.7) channel. 

Electrons are required to have at least 90% of their en- 
ergy deposited in the electromagnetic (EM) calorimeter 
and a distribution for EM shower consistent with that 
expected for an electron. They are further required to be 
isolated in both the calorimeter and the tracker. A neu- 
ral network (NN) multivariate discriminant [10j, formed 
from the parameters of the EM shower and the track as- 
sociated with the electron candidate, as well as central 
preshower detector information, is used to discriminate 
electrons from jets. For electrons with pt = 40 GeV, the 
identification efficiency is s=s 82%. 

Muons are identified as track segments in the muon 
detector that match tracks found in the tracking system. 
They must be synchronous with the beam crossing time 
to reject background from cosmic rays. Muons are also 
required to be isolated in both the calorimeter and the 
tracker. The identification efficiency for muons with pt = 
40 GeV is w 79%. 

Photons are identified in the central calorimeter (CC) 
and are required to be separated from leptons and jets by 
AR = ^(Ar/) 2 + (A0) 2 > 0.7 |ll|]. Additional require- 
ments are applied on the fraction of energy deposited in 
the EM calorimeter and on isolation in both the calorime- 
ter and the tracker. The shower width in the third layer 
(EM3) must be consistent with that of a photon. To sup- 
press electrons misidentified as photons, the candidates 
must not be spatially matched to a track or to energy de- 
positions in the silicon microstrip or central fiber trackers 
that lie along the trajectory of an electron [12j. Fur- 
ther rejection of jets is achieved with a NN discriminant 
similar to that used for electron selection. The average 
identification efficiency for photons with pt = 40 GeV is 
« 75%. 

The $t is the negative of the vectorial sum of trans- 
verse components of energy depositions in the calorime- 
ter, corrected for identified photons, electrons, and 
muons. Jet energies are calibrated using transverse en- 
ergy balance in photon+jet and dijet events [13|, and 
these corrections are propagated to the calculation of $ T . 

To select Z7 + $ T events, we first require at least 
two leptons that are consistent with originating from 
Z — > e + e~ or /i + \x~~ decay. Each lepton must have 
Pt > 15 GeV, with one electron (muon) having pr > 
25 (20) GeV. The two leptons must have opposite charge 
and an invariant mass Mm within the Z-mass windows 
of 78-104 GeV and 65-115 GeV for the ee and \x\x chan- 
nels, respectively. A total number of 261,964 (306,541) 
ee (p-fJ-) candidates satisfy these criteria. We require 
at least one isolated photon with p T > 30 GeV in the 
event. To reduce background from photons radiated by 
the two leptons, we require a three-body invariant mass 



M(££j) > 120 GeV, which results in a total number 
of 78 (91) ee7 (fifJfy) candidates. The GMSB signal is 
expected in the region of large $ T . We therefore re- 
quire $t > 30 (40) GeV in the electron (muon) channel. 
To remove events with spurious $ T due to poorly re- 
constructed muons, we require that A<j>{$ T ,ni) < 2.85 
where fj,% is the highest-pr muon. The f] T significance, 
a likelihood discriminant based on the ratio of $? an d 
its uncertainty, is required to be > 5. These selections 
optimize sensitivity to signal. No data is selected in the 
ee7 + $ T final state, and a single event is selected in the 
fifi'y + $t final state. 

The background to the Z7 + $ T signal arises from 
instrumental backgrounds caused by mismeasured flxi 
misidentified leptons or misidentified jets in Z7, Z+jets, 
WW, WZ, ZZ, W + X,tt and multijet processes. The 
backgrounds are either estimated using control samples 
in data or using Monte Carlo (MC) simulated events pro- 
cessed using a detailed GEANT-based simulation 14| of 
the DO detector response and overlaid with data from 
random beam crossings. The simulation is corrected for 
lepton identification efficiencies and energy resolutions 
observed in data. 

The SM Z7 process is the dominant source of back- 
ground. It is estimated using PYTHIA [15|. The photon 
Pt spectrum from PYTHIA for initial state radiation (ISR) 
is corrected for QCD and electroweak next-to-leading 
order (NLO) effects using the MC event generator of 
Ref. (lfjj. The contribution from final state radiation in 
data is determined by fitting the M{££) distribution of 
Z7 MC events to data in the range p T > 10 GeV and 
$ T < 30 GeV and is found to be very small because of 
the requirements on AR(£,^), p T , M(££j), and $ T . We 
estimate the Z7 contribution in the signal region to be 
0.23 ± 0.05 (stat) and 0.43 ± 0.05 (stat) events in the 
electron and muon channels, respectively. 

Background from Z+jets events can enter the sam- 
ple if a jet is misidentified as a photon and $ T is large. 
Two data-driven methods are used to estimate this back- 
ground. In the first method, we select an orthogonal sam- 
ple of events with at least two electrons or two muons and 
with a jet passing all photon acceptance criteria except 
failing either the requirements on tracker isolation or on 
shower width in EM3. The Z+jets background is then 
estimated by scaling this sample by an 7/— dependent fac- 
tor /. This factor / is the ratio of the probability for 
a jet to satisfy full photon-identification criteria to the 
probability to fail tracker isolation or shower width re- 
quirements. It is measured using dijet data as a func- 
tion of 77 and Et, yielding typical values of 0.08 to 0.16 
with uncertainties of 10%. In the second method, the 
Z+jets background is estimated by fitting the sum of 
the NN templates for photons and photon-like jets to the 
observed photon NN distribution. Templates of the NN 
distributions are obtained from simulations of photons 
and separately of jets, as the NN for data is found to be 



TABLE I. Cross sections <r p for the production of pairs of 
lightest neutralinos Xi vla cascade decay, branching fractions 
of x? to 7G (B 7 ) and to ZG (Bz), and the lightest neutralino 
mass M-o used in this analysis, which is parametrized by 

the breaking scale A. The \i a ^ so decays to Higgs+G and 
to nonresonant £ + £~G, which dominate the remaining decays 
for large and small A, respectively. Also given are the ob- 
served (expected) 95% C.L. upper limits on the production 
cross section. 
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well modeled by MC [lOj. The results from these two 
methods are consistent within their statistical uncertain- 
ties, and the first method is used since it yields smaller 
uncertainties. The resulting estimates of the Z+jets con- 
tribution in the signal region are 0.09 ± 0.08 (stat) and 
0.17 ± 0.16 (stat) in the electron and muon channels, re- 
spectively. 

The multijet contribution to the background for Z — > 
II candidates is estimated by fitting the M« distribu- 
tion using templates from jet-rich data and MC simu- 
lated Z — >• It events. Using weighted jet- rich data, the 
contribution in the signal region is found to be negligible. 

The SM backgrounds from WW, WZ, ZZ, and ti pro- 
duction are estimated using MC simulations. The j5 T can 
be substantial in such events, but none of these back- 
grounds are sources of isolated, high-pj photons. The 
contribution from ti events is minimized by the require- 
ments on Ma. 

The GMSB signal is modeled with the PYTHIA leading- 
order (LO) MC event generator using supersymmetric 
particle spectra calculated in ISAJET 17J. The A pa- 



rameter is varied from 70 TeV to 95 TeV, in steps of 
5 TeV, and used to compute an MSSM particle mass 
spectrum and a set of branching ratios. The LO signal 
cross sections are scaled to match the NLO prediction 
from PROSPINO [18j. The inclusive cross section for the 
pair production of Xi from cascade decays is 618 fb for 
A = 70 TeV and decreases to 106 fb for A = 95 TeV. 
The fraction of Xi — > ZG decays (&z) increases with A, 
reaching 50% at A sa 85 TeV. Cross sections and branch- 
ing fractions are given in Table Q] At larger A values, ZG 
is the main decay mode for Xi- For the full event selec- 
tion, the overall product of acceptance and efficiency is 
7.7 (5.1)% at A = 70 TeV and increases to 11.2 (8.6)% 
for A = 95 TeV in the electron (muon) channel. 

The expected signal yield for A = 80 and 90 TeV and 



the estimated SM backgrounds are summarized in Ta- 
ble [TTJ The total background is expected to be 0.5 ± 0.1 
and 0.7 ± 0.4 events in the eej + f T and /z/ry + $? 
channels, respectively. The number of observed events 
is consistent with these expectations. The comparison 
between data and SM MC predictions for the f T distri- 
butions after selecting Z7 events is given in Fig. [1] along 
with the signal expectation. Good agreement between 
data and SM background is observed for both eej and 
/x/X7 channels. 

The systematic uncertainties that affect the signal and 
SM backgrounds include theoretical and experimental 
sources. The uncertainties on the theoretical cross sec- 
tion for diboson and ti processes are 6% and 10%, re- 
spectively. The uncertainty on the measured luminosity 
is 6.1% [9j and is applied to the SM background estima- 
tions based on MC simulation. The uncertainty on elec- 
tron identification efficiency is 1% in the CC region and 
increases to 4% in the end-cap calorimeter. The system- 
atic uncertainties on muon identification include 1.0% for 
reconstruction, 1.1% for tracking efficiency, and 0.5% for 
isolation. The photon identification uncertainty is 2.7%. 
The uncertainties from the jet energy scale are estimated 
to be 1% for signal and 4% for the backgrounds. The 
uncertainty on the momentum resolution for muons is 
reflected in an uncertainty of « 100% in the signal re- 
gion $ T > 40 GeV on the estimate of the background 
from Z(fifj,) + 7. 

To improve the sensitivity for X? detection at the cost 
of a stronger dependence on the specifics of the GMSB 
model, we also use a BDT multivariate technique to dis- 
criminate between SM background and signal [19J. The 
output is a discriminant that is shifted toward +1 for sig- 
nal, and strongly peaked near —1 for background events. 

The BDT is trained on a randomly selected collection 
of signal and background MC events, Z+jet background 
candidates from data, and a signal assuming A = 90 TeV. 
The training samples require a leading lepton of px > 
25 (20) GeV, a second lepton of p T > 15 GeV, pj > 20 
GeV, M(£l-y) > 120 GeV, $ T > 15 GeV and M(U) > 
70 (65) GeV in the electron (muon) channel. A set of 
14 sensitive variables, well modeled by the simulation, 
is used to form the BDT discriminant. The variables 
include transverse momenta of the two leptons, photon, 
dilepton system, and dilepton+photon system, as well as 
$ T and M(££j). The expected signal and background 
yields are estimated from events independent of the set 
used for training. The data is found consistent with the 
SM background prediction as seen in Fig. [5] and Table |TT] 
(for BDT > 0.8), and no evidence is observed for a GMSB 
neutralino NLSP. 

Limits on the production cross section of X1X1 using 
the benchmark model are derived using a Poisson log- 
likelihood ratio as test statistic, combining results from 
the electron and muon channels. Pseudo-experiments 
are generated according to the background-only and sig- 



TABLE II. Number of observed and expected events for the restrictive criteria defining the signal region and for less stringent 
requirements that are followed by a selection on BDT output defining an alternative signal region. The first uncertainty is 
statistical and the second is systematic. The contributions from Z+jets for the BDT-analyses are found to be negligible. 
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3.28 ±0.09 ±0.24 
1.48 ±0.03 ±0.11 
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1.73 ±0.05 ±0.21 
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ZZ 



0.23 ±0.05 ±0.02 
0.09 ±0.08 ±0.01 
0.13±0.05±0.01 
0.05 ±0.01 ±0.01 
0.50 ±0.11 ±0.03 
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0.14 ±0.03 ±0.02 
0.43 ±0.12 ±0.03 



0.43 ±0.05 ±0.40 
0.17±0.16±0.02 
0.08 ±0.03 ±0.01 
0.04 ±0.01 ±0.01 
0.71±0.17±0.40 



0.10 ±0.03 ±0.20 

0.16±0.19±0.02 
0.05 ±0.02 ±0.01 
0.31 ±0.10 ±0.20 
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FIG. 1. Distribution of 1JZ T for Z^j events in the (a) eey channel and (b) ^7 channel before requiring $ T significance > 5. 
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FIG. 2. Distribution of BDT output in the ee'y channel and 
fi/j,y channel for background only, background with a A = 
85 TeV signal added, and for data. The total background 
uncertainties are indicated as shaded bands. 



FIG. 3. Limit on the cross section for Zy + fi T production 
as a function of A (lower horizontal axis) and M(xl) (upper 
horizontal axis) at 95% C.L. combined for the ee'y and fifiy 
channels. The NLO cross section from theory is overlaid. 



nal±background hypotheses, and systematic uncertain- 
ties are accounted for by integrating over uncertainties 
parametrized as Gaussian. The limits on cross sec- 
tions are evaluated using the modified frequentist ap- 



proach [20(. Data and background estimates are studied 
in four bins of BDT output, and values from the most 
signal-like bin (BDT > 0.8) are shown in Table HH The 
95% C.L. upper limit on the cross section using the BDT 



discriminant is shown in Fig. [3l together with the ex- 
pected limit and the 1 and 2 standard deviation (SD) 
uncertainty bands. The 95% C.L. limits on a p are also 
given in Table Q] Scales of A < 87 TeV are excluded at 
95% C.L. which corresponds to M(xT) < 151 GeV. 

In summary, we present the first search for a SUSY 
signature in events containing Z7 + $ T final states using 
6.2 fb _1 of integrated luminosity collected by the DO ex- 
periment inpp collisions at y^ =1.96 TeV. The signature 
corresponds to a GMSB model where pairs of neutralino 
NLSPs are either produced promptly or from decays of 
other supersymmetric particles in pp collisions and then 
decay to either ZG or jG. In the expected signal re- 
gion we observe no event in the eej + $t an( i one event 
in the /i/17 + $ T channels, where the SM background is 
expected to be 1.21 ±0.45 combined. Employing a multi- 
variate selection process and combining the results from 
both channels, the specific neutralino NLSP model is ex- 
cluded at the 95% C.L. for A < 87 TeV, corresponding 
to neutralino masses of M(x°) < 151 GeV. 
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